Strategy proposal using QSAR models to approach mutagenicity assessment of non intentionally added substances in recycled plastic resins.

CONTEXT: Transition to the use of recycled plastics raises an issue concerning safety assessment of Non Intentionally Added Substances (NIAS). To assess the mutagenic potential of the recycled polyethylene impurities and to evaluate the need to perform in vitro assays on recycled resins, this study lies in identifying existing NIAS associated with recycled Low/High Density Polyethylene and assessing the mutagenicity data-gaps by employing in silico tools. METHODS: Quantitative Structure-Activity Relationship (QSAR) models predicting Ames mutagenicity were selected from literature, then NIAS were run to 1/evaluate performances of each model, 2/apply a QSAR strategy on the NIAS molecular space and address data-gaps. RESULTS: Among the 165 NIAS identified, experimental Ames results were not found for 50 substances while the substances with experimental data were predominantly negatives. No individual model was able to predict all NIAS due to applicability domain limitations. Taking into account 1/calculated performances, 2/availability of applicability domain, 3/description of the Training Set, an Integrated Strategy was founded including Sarpy, Consensus and Protox to extend the applicability domain. CONCLUSION & PERSPECTIVES: Existing data and predictions generated by this strategy suggest a low mutagenic potential of NIAS. Further investigation is needed to explore other genotoxicity mechanisms.

Assessment of the performance of the Ames MPF™ assay: A multicenter collaborative study with six coded chemicals.

The Ames MPF™ is a miniaturized, microplate fluctuation format of the Ames test. It is a standardized, commercially available product which can be used to assess mutagenicity in Salmonella and E. coli strains in 384-well plates using a color change-based readout. Several peer-reviewed comparisons of the Ames MPF™ to the Ames test in Petri dishes confirmed its suitability to evaluate the mutagenic potential of a variety of test items. An international multicenter study involving seven laboratories tested six coded chemicals with this assay using five bacterial strains, as recommended by the OECD test guideline 471. The data generated by the participating laboratories was in excellent agreement (93%), and the similarity of their dose response curves, as analyzed with sophisticated statistical approaches further confirmed the suitability of the Ames MPF™ assay as an alternative to the Ames test on agar plates, but with advantages with respect to significantly reduced amount of test substance and S9 requirements, speed, hands-on time and, potentially automation.

Mutagenicity assessment of high-power 1.6-THz pulse laser radiation.

The effect of terahertz (THz) radiation has been studied in medicine. However, there is a lack of scientific information regarding its possible mutagenicity. Therefore, the present study aimed to assess the mutagenicity of 1.6 THz laser irradiation. The Ames test was conducted using five bacterial tester strains. The bacteria were subjected to (i) 1.6 THz laser irradiation at 3.8 mW/cm2 for 60 min using a tabletop THz pulse laser system, (ii) ultraviolet irradiation, (iii) treatment with positive control chemicals (positive control) or (iv) treatment with the solvent used in the positive control (negative control). After treatment, the bacterial suspensions were cultured on minimal glucose agar to determine the number of revertant colonies. In addition, the comet assay was performed using fibroblasts (V79) to assess possible DNA damage caused by the THz laser irradiation. The Ames test demonstrated that the THz laser irradiation did not increase the number of revertant colonies compared to that in the negative control group, whereas the ultraviolet irradiation and positive control treatment increased the number of revertant colonies. Thus, 1.6 THz laser irradiation is unlikely to be mutagenic. The comet assay additionally suggests that the THz laser irradiation unlikely induce cellular DNA damage.

Bacterial worth in genotoxicity assessment studies.

Bacterial-based genotoxicity test systems play a significant role in the detection and evaluation of genotoxicity in vitro and have gained importance due to attributes like wide applicability, speed, high sensitivity, good reproducibility, and simplicity. The Salmonella microsomal mutagenicity assay was created by Ames and colleagues at the beginning of the 1970s, and it was based on the fundamental notion that in auxotrophic bacterial strains with inhibited growth, a mutant gene would revert to its original state on exposure to genotoxicants. This is the most successful and widely used in vitro genotoxicity test. Later, a number of additional test systems that incorporated DNA repair mechanisms including the bacterial SOS response were created. Genetic engineering has further provided significant advancement in these test systems with the development of highly sophisticated bacterial tester strains with significantly increased sensitivity to evaluate the chemical nature of hazardous substances and pollutants. These bacterial bioassays render an opportunity to detect the defined effects of compounds at the molecular level. In this review, all the aspects related to the bacterial system in genotoxicity assessment have been summarized and their role is elaborated concerning real-time requirements and future perspectives.

EnzyHTP Computational Directed Evolution with Adaptive Resource Allocation.

Directed evolution facilitates enzyme engineering via iterative rounds of mutagenesis. Despite the wide applications of high-throughput screening, building "smart libraries" to effectively identify beneficial variants remains a major challenge in the community. Here, we developed a new computational directed evolution protocol based on EnzyHTP, a software that we have previously reported to automate enzyme modeling. To enhance the throughput efficiency, we implemented an adaptive resource allocation strategy that dynamically allocates different types of computing resources (e.g., GPU/CPU) based on the specific need of an enzyme modeling subtask in the workflow. We implemented the strategy as a Python library and tested the library using fluoroacetate dehalogenase as a model enzyme. The results show that compared to fixed resource allocation where both CPU and GPU are on-call for use during the entire workflow, applying adaptive resource allocation can save 87% CPU hours and 14% GPU hours. Furthermore, we constructed a computational directed evolution protocol under the framework of adaptive resource allocation. The workflow was tested against two rounds of mutational screening in the directed evolution experiments of Kemp eliminase (KE07) with a total of 184 mutants. Using folding stability and electrostatic stabilization energy as computational readout, we identified all four experimentally observed target variants. Enabled by the workflow, the entire computation task (i.e., 18.4 µs MD and 18,400 QM single-point calculations) completes in 3 days of wall-clock time using ∼30 GPUs and ∼1000 CPUs.

Assessment of DNA mutagenicity induced by He-Ne laser using Salmonella typhimurium strains.

Helium-neon (He-Ne) laser mutagenesis is widely used in microbiology and plant breeding. In this study, two frameshift mutant representative strains of Salmonella typhimurium TA97a and TA98 and two base pair substitution types TA100 and TA102 were employed as model microorganisms to assess DNA mutagenicity induced by He-Ne laser (3 J·cm-2·s-1, 632.8 nm) for 10, 20, and 30 min. The results revealed that the optimal laser application was 6 h in the mid-logarithmic growth stage. Low-power He-Ne laser for short treatment inhibited cell growth, and continued treatment stimulated the metabolism. The effects of the laser on TA98 and TA100 were the most prominent. Sequencing results from 1500 TA98 revertants showed that there were 88 insertion and deletion (InDel) types in the hisD3052 gene, of which the InDels unique to laser were 21 more than that of the control. Sequencing results from 760 TA100 revertants indicated that laser treatment created Pro (CCC) in the product of the hisG46 gene more likely to be replaced by His (CAC) or Ser (TCC) than by Leu (CTC). Two unique non-classical base substitutions, CCC → TAC and CCC → CAA, also appeared in the laser group. These findings will provide a theoretical basis for further exploration of laser mutagenesis breeding. KEY POINTS: • Salmonella typhimurium served as model organism for laser mutagenesis study. • Laser promoted the occurrence of InDels in the hisD3052 gene of TA98. • Laser promoted the occurrence of base substitution in the hisG46 gene of TA100.

Enabling low-cost and robust essentiality studies with high-throughput transposon mutagenesis (HTTM).

Transposon-insertion sequencing (TIS) methods couple high density transposon mutagenesis with next-generation sequencing and are commonly used to identify essential or important genes in bacteria. However, this approach can be work-intensive and sometimes expensive depending on the selected protocol. The difficulty to process a high number of samples in parallel using standard TIS protocols often restricts the number of replicates that can be performed and limits the deployment of this technique to large-scale projects studying gene essentiality in various strains or growth conditions. Here, we report the development of a robust and inexpensive High-Throughput Transposon Mutagenesis (HTTM) protocol and validate the method using Escherichia coli strain BW25113, the parental strain of the KEIO collection. HTTM reliably provides high insertion densities with an average of one transposon every ≤20bp along with impressive reproducibility (Spearman correlation coefficients >0.94). A detailed protocol is available at protocol.io and a graphical version is also included with this article.

Mutation resource of Samba Mahsuri revealed the presence of high extent of variations among key traits for rice improvement.

To create novel variants for morphological, physiological, and biotic stress tolerance traits, induced mutations were created using Ethyl Methane Sulphonate (EMS) in the background of Samba Mahsuri (BPT 5204), a popular and mega rice variety of India. A population derived from 10, 500 M1 plants and their descendants were phenotyped for a wide range of traits leading to the identification of 124 mutants having variations in key agro-morphological traits, and 106 mutants exhibiting variation for physiological traits. Higher yield is the ultimate goal of crop improvement and we identified 574 mutants having higher yield compared to wild type by having better yield attributing traits. Further, a total of 50 mutants showed better panicle exertion phenotypes as compared to Samba Mahsuri leading to enhancement of yield. Upon rigorous screening for three major biotic stresses, 8 mutants showed enhanced tolerance for yellow stem borer (YSB), and 13 different mutants each showed enhanced tolerance for sheath blight (ShB) and bacterial leaf blight (BLB), respectively. In addition, screening at multiple locations that have diverse field isolates identified 3, 3, and 5 lines for tolerance to ShB, YSB and BLB, respectively. On the whole, 1231 desired mutant lines identified at M2 were forwarded to an advanced generation (M5). PCR based allele mining indicated that the BLB tolerant mutants have a different allele than the reported alleles for well-known genes affecting bacterial blight resistance. Whole genome re-sequencing revealed substantial variation in comparison to Samba Mahsuri. The lines showing enhanced tolerance to important biotic stresses (YSB, ShB and BLB) as well as several economically important traits are unique genetic resources which can be utilized for the identification of novel genes/alleles for different traits. The lines which have better agronomic features can be used as pre-breeding lines. The entire mutant population is maintained as a national resource for genetic improvement of the rice crop.

The codification of hazard and its impact on the hazard versus risk controversy.

The long running controversy about the relative merits of hazard-based versus risk-based approaches has been investigated. There are three levels of hazard codification: level 1 divides chemicals into dichotomous bands of hazardous and non-hazardous; level 2 divides chemicals into bands of hazard based on severity and/or potency; and level 3 places each chemical on a continuum of hazard based on severity and/or potency. Any system which imposes compartments onto a continuum will give rise to issues at the boundaries, especially with only two compartments. Level 1 schemes are only justifiable if there is no variation in severity, or potency or if there is no threshold. This is the assumption implicit in GHS/EU classification for carcinogenicity, reproductive toxicity and mutagenicity. However, this assumption has been challenged. Codification level 2 hazard assessments offer a range of choices and reduce the built-in conflict inherent in the level 1 process. Level 3 assessments allow a full range of choices between the extremes and reduce the built-in conflict even more. The underlying reason for the controversy between hazard and risk is the use of level 1 hazard codification schemes in situations where there are ranges of severity and potency which require the use of level 2 or level 3 hazard codification. There is not a major difference between level 2 and level 3 codification, and they can both be used to select appropriate risk management options. Existing level 1 codification schemes should be reviewed and developed into level 2 schemes where appropriate.

Assessment of mutations induced by cold atmospheric plasma jet treatment relative to known mutagens in Escherichia coli.

The main bactericidal components of cold atmospheric plasma (CAP) are thought to be reactive oxygen and nitrogen species (RONS) and UV-radiation, both of which have the capacity to cause DNA damage and mutations. Here, the mutagenic effects of CAP on Escherichia coli were assessed in comparison to X- and UV-irradiation. DNA damage and mutagenesis were screened for using a diffusion-based DNA fragmentation assay and modified Ames test, respectively. Mutant colonies obtained from the latter were quantitated and sequenced. CAP was found to elicit a similar mutation spectrum to X-irradiation, which did not resemble that for UV implying that CAP-produced RONS are more likely the mutagenic component of CAP. CAP treatment was also shown to promote resistance to the antibiotic ciprofloxacin. Our data suggest that CAP treatment has mutagenic effects that may have important phenotypic consequences.

In vitro mutagenicity assessment of fried meat-based food from mass catering companies.

The current article aimed to evaluate the in vitro mutagenicity of ten fried meat-based food extracts obtained from different catering companies from Navarra (Spain). A miniaturized 6-well version of the Ames test in Salmonella typhimurium TA98, and the in vitro micronucleus test (OECD TG 487) in TK6 cells were performed. None of the ten extracts of fried meat-based food induced gene mutations in S. typhimurium TA98 with or without metabolic activation, but five induced chromosomal aberrations after 24 h treatment of TK6 without metabolic activation. More studies are needed to check the biological relevance of these in vitro studies.

Applied mutagenesis could improve economically important traits in bitter gourd (Momordica charantia L.).

Mutants with unique characters have played a key role in discovery of gene, mapping, functional genomics and breeding in many vegetable crops, but information on bitter gourd is lacking. Induction of mutation by gamma rays (Co60 source) at five different doses (50 Gy, 100 Gy, 150 Gy, 200 Gy and 250 Gy) was studied in four widely divergent bitter gourd genotypes BG-1346501, Meghna-2, Special Boulder and Selection-1 in M1 generation. Reduction in seed germination percentage, vine length and pollen fertility occurred in M1 generation with the increasing doses of mutagens. LD50 dose for BG-1346501, Meghna-2, Special Boulder and Selection-1 corresponded to 290.76 Gy, 206.12 Gy, 212.81 Gy and 213.49 Gy ᵞ radiation, respectively suggested low to medium doses (200-250 Gy) of gamma rays would be helpful in producing useful and exploitable mutants for further breeding. No remarkable effect of ᵞ radiation on fruit physicochemical characters in M1 generation were observed. M2 generation, raised from two widely divergent genotypes, BG-1346501 and Meghna-2, were screened critically and observed no significant reduction in seed germination and pollen viability, however little damage occurred particularly in vine length. There is possibility of isolating segregates in M2 generation with enhanced nutrient contents at low radiation dose. Highest mutation frequency resulted by treating Meghna-2 at 200 Gy and BG-1346501 at 100 Gy. Both genotype and mutagenic doses influenced mutagenic effectiveness. Spectrum of mutation was very low; number of putative mutants isolated from M2 generation was five in Meghna-2 and three in BG-1346501. Among six putative macro mutants isolated from M3 generation, we could identify two putative mutants, namely Meghna-2 with gynoecious sex form and BG-1346501 with high charantin, appreciable ß-carotene and high ascorbic acid contents having ample promise for further utilization in bitter gourd breeding after critical testing in subsequent generations for estimation of genetic gain and trait heritability to confirm the mutant stability.

Improving sequence-based modeling of protein families using secondary-structure quality assessment.

MOTIVATION: Modeling of protein family sequence distribution from homologous sequence data recently received considerable attention, in particular for structure and function predictions, as well as for protein design. In particular, direct coupling analysis, a method to infer effective pairwise interactions between residues, was shown to capture important structural constraints and to successfully generate functional protein sequences. Building on this and other graphical models, we introduce a new framework to assess the quality of the secondary structures of the generated sequences with respect to reference structures for the family. RESULTS: We introduce two scoring functions characterizing the likeliness of the secondary structure of a protein sequence to match a reference structure, called Dot Product and Pattern Matching. We test these scores on published experimental protein mutagenesis and design dataset, and show improvement in the detection of nonfunctional sequences. We also show that use of these scores help rejecting nonfunctional sequences generated by graphical models (Restricted Boltzmann Machines) learned from homologous sequence alignments. AVAILABILITY AND IMPLEMENTATION: Data and code available at https://github.com/CyrilMa/ssqa. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Massively parallel assessment of human variants with base editor screens.

Understanding the functional consequences of single-nucleotide variants is critical to uncovering the genetic underpinnings of diseases, but technologies to characterize variants are limiting. Here, we leverage CRISPR-Cas9 cytosine base editors in pooled screens to scalably assay variants at endogenous loci in mammalian cells. We benchmark the performance of base editors in positive and negative selection screens, identifying known loss-of-function mutations in BRCA1 and BRCA2 with high precision. To demonstrate the utility of base editor screens to probe small molecule-protein interactions, we screen against BH3 mimetics and PARP inhibitors, identifying point mutations that confer drug sensitivity or resistance. We also create a library of single guide RNAs (sgRNAs) predicted to generate 52,034 ClinVar variants in 3,584 genes and conduct screens in the presence of cellular stressors, identifying loss-of-function variants in numerous DNA damage repair genes. We anticipate that this screening approach will be broadly useful to readily and scalably functionalize genetic variants.

Assessment of efficacy of mutagenesis of gamma-irradiation in plant height and days to maturity through expression analysis in rice.

Introduction of semi-dwarfism and early maturity in rice cultivars is important to achieve improved plant architecture, lodging resistance and high yield. Gamma rays induced mutations are routinely used to achieve these traits. We report the development of a semi-dwarf, early maturing and high-yielding mutant of rice cultivar 'Improved White Ponni', a popular cosmopolitan variety in south India preferred for its superior grain quality traits. Through gamma rays induced mutagenesis, several mutants were developed and subjected to selection up to six generations (M6) until the superior mutants were stabilized. In the M6 generation, significant reduction in days to flowering (up to 11.81% reduction) and plant height (up to 40% reduction) combined with an increase in single plant yield (up to 45.73% increase) was observed in the mutant population. The cooking quality traits viz., linear elongation ratio, breadthwise expansion ratio, gel consistency and gelatinization temperature of the mutants were similar to the parent variety Improved White Ponni. The genetic characterization with SSR markers showed variability between the semi-dwarf-early mutants and the Improved White Ponni. Gibberellin responsiveness study and quantitative real-time PCR showed a faulty gibberellin pathway and epistatic control between the genes such as OsKOL4 and OsBRD2 causing semi-dwarfism in a mutant. These mutants have potential as new rice varieties and can be used as new sources of semi-dwarfism and earliness for improving high grain quality rice varieties.

Assessment of Adenosine Triphosphatase Phospholipid Transporting 8B1 (ATP8B1) Function in Patients With Cholestasis With ATP8B1 Deficiency by Using Peripheral Blood Monocyte-Derived Macrophages.

Adenosine triphosphatase phospholipid transporting 8B1 (ATP8B1) deficiency, an ultrarare autosomal recessive liver disease, includes severe and mild clinical forms, referred to as progressive familial intrahepatic cholestasis type 1 (PFIC1) and benign recurrent intrahepatic cholestasis type 1 (BRIC1), respectively. There is currently no practical method for determining PFIC1 or BRIC1 at an early disease course phase. Herein, we assessed the feasibility of developing a diagnostic method for PFIC1 and BRIC1. A nationwide Japanese survey conducted since 2015 identified 25 patients with cholestasis with ATP8B1 mutations, 15 of whom agreed to participate in the study. Patients were divided for analysis into PFIC1 (n = 10) or BRIC1 (n = 5) based on their disease course. An in vitro mutagenesis assay to evaluate pathogenicity of ATP8B1 mutations suggested that residual ATP8B1 function in the patients could be used to identify clinical course. To assess their ATP8B1 function more simply, human peripheral blood monocyte-derived macrophages (HMDMs) were prepared from each patient and elicited into a subset of alternatively activated macrophages (M2c) by interleukin-10 (IL-10). This was based on our previous finding that ATP8B1 contributes to polarization of HMDMs into M2c. Flow cytometric analysis showed that expression of M2c-related surface markers cluster of differentiation (CD)14 and CD163 were 2.3-fold and 2.1-fold lower (95% confidence interval, 2.0-2.5 for CD14 and 1.7-2.4 for CD163), respectively, in patients with IL-10-treated HMDMs from PFIC1 compared with BRIC1. Conclusion: CD14 and CD163 expression levels in IL-10-treated HMDMs may facilitate diagnosis of PFIC1 or BRIC1 in patients with ATP8B1 deficiency.

A forward selection algorithm to identify mutually exclusive alterations in cancer studies.

Mutual exclusivity analyses provide an effective tool to identify driver genes from passenger genes for cancer studies. Various algorithms have been developed for the detection of mutual exclusivity, but controlling false positive and improving accuracy remain challenging. We propose a forward selection algorithm for identification of mutually exclusive gene sets (FSME) in this paper. The method includes an initial search of seed pair of mutually exclusive (ME) genes and subsequently including more genes into the current ME set. Simulations demonstrated that, compared to recently published approaches (i.e., CoMEt, WExT, and MEGSA), FSME could provide higher precision or recall rate to identify ME gene sets, and had superior control of false positive rates. With application to TCGA real data sets for AML, BRCA, and GBM, we confirmed that FSME can be utilized to discover cancer driver genes.

Filling the gap between risk assessment and molecular determinants of tumor onset.

In the past two decades, a ponderous epidemiological literature has causally linked tumor onset to environmental exposure to carcinogens. As consequence, risk assessment studies have been carried out with the aim to identify both predictive models of estimating cancer risks within exposed populations and establishing rules for minimizing hazard when handling carcinogenic compounds. The central assumption of these works is that neoplastic transformation is directly related to the mutational burden of the cell without providing further mechanistic clues to explain increased cancer onset after carcinogen exposure. Nevertheless, in the last few years, a growing number of studies have implemented the traditional models of cancer etiology, proposing that neoplastic transformation is a complex process in which several parameters and crosstalk between tumor and microenvironmental cells must be taken into account and integrated with mutagenesis. In this conceptual framework, the current strategies of risk assessment that are solely based on the 'mutator model' require an urgent update and revision to keep pace with advances in our understanding of cancer biology. We will approach this topic revising the most recent theories on the biological mechanisms involved in tumor formation in order to envision a roadmap leading to a future regulatory framework for a new, protective policy of risk assessment.

Comprehensive assessment of TP53 loss of function using multiple combinatorial mutagenesis libraries.

The diagnosis of somatic and germline TP53 mutations in human tumors or in individuals prone to various types of cancer has now reached the clinic. To increase the accuracy of the prediction of TP53 variant pathogenicity, we gathered functional data from three independent large-scale saturation mutagenesis screening studies with experimental data for more than 10,000 TP53 variants performed in different settings (yeast or mammalian) and with different readouts (transcription, growth arrest or apoptosis). Correlation analysis and multidimensional scaling showed excellent agreement between all these variables. Furthermore, we found that some missense mutations localized in TP53 exons led to impaired TP53 splicing as shown by an analysis of the TP53 expression data from the cancer genome atlas. With the increasing availability of genomic, transcriptomic and proteomic data, it is essential to employ both protein and RNA prediction to accurately define variant pathogenicity.